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TRANSFORMATION OF APPARENT OCEAN WAVE SPECTRA 
OBSERVED FROM AN AIRCRAFT SENSOR PLATFORM 


Lamont R. Poole 
Langley Research Center 

SUMMARY 

The problem considered herein was the transformation of a unidirectional apparent 
ocean wave spectrum observed from an aircraft sensor platform into the true spectrum, 
or the spectrum which would result from observation of the same wave field from a sta- 
tionary platform. Spectral transformation equations valid for all depth zones were devel- 
oped in terms of the linear wave dispersion relationship, and an iterative solution for a 
specified average water depth was outlined. Results obtained by using this iterative solu- 
tion indicated that different transformed spectra can be obtained by varying the average 
water depth in the transformation of a single apparent spectrum. The differences were 
most significant when much of the energy density was expected, a priori, to exist at rela- 
tively low true frequencies. Varying the average depth led to a redistribution of energy 
density among the various frequency bands of the transformed spectrum while the total 
energy density in the spectrum was preserved. The necessity of seeking an iterative 
solution to the transformation equations was emphasized in an appendix in which errors 
resulting from the use of certain approximations were discussed. 

INTRODUCTION 

With man^s increasing usage of the continental shelf regions of the world, the need 
for techniques to monitor and study broad-scale ocean phenomena such as surface waves 
and currents has been emphasized. One tool which could potentially play an important 
role in an overall coastal-zone management program is the wave refraction model, such 
as that of reference 1. A wave refraction model, in theory, predicts the behavior of ocean 
surface waves as they cross the continental shelf and impinge on the shoreline. If the 
model is of a sufficiently large spatial scale that regional differences in wave environment 
can be highlighted, the process of site selection for offshore and coastal activities is 
aided. 

One must consider, however, the fact that the wave refraction model is a mathemat- 
ical simplification of a complex physical process. In this light, experimental verification 
of model results is imperative before the wave refraction model can be used in an opera- 



tional sense for prediction and decision making. The broad-scale experimental data on 
the coastal wave environment required for model verification can best be obtained through 
the use of aircraft carrying instruments such as a laser profilometer or radar altimeter 
to measure small-scale variations in sea-surface elevation. (See ref. 2, for example.) 

An important computed product of such surface elevation data is the wave energy density 
(energy per unit area) spectrum, which in its most common form displays the wave energy 
density as a function of wave frequency. Such energy density spectra could be obtained 
for a number of deep-water and nearshore locations along the flight path of the aircraft. 

The accuracy of the wave refraction model could then be assessed by using the model, as 
in reference 3, to simulate the variations in wave energy density in discrete bands of the 
experimentally derived spectra as the wave field propagates from deep-water to nearshore 
locations. 

A problem which must be considered with the use of such remotely sensed surface 
elevation data is that the observations have been made from a rapidly moving reference 
frame, the aircraft. A frequency spectrum of wave energy computed from these data 
accordingly is valid in the moving reference frame but provides a distorted estimate of 
the true frequency spectrum in a fixed reference frame. The problem of mapping this 
apparent frequency spectrum into its counterpart in a fixed reference frame was first dis- 
cussed in a classic paper by St. Denis and Pierson (ref. 4). The basic relationships given 
in reference 4 were discussed from a more practical computational viewpoint by Cartwright 
(ref. 5) in a study of the output of a wave recorder on board a moving ship. More recently 
these relationships have been used in a study of wave generation and growth (ref. 6) per- 
formed by using airborne instrumentation in a mode similar to that desired for a refrac- 
tion model verification study. 

In the previously mentioned references concerning transformation of apparent spec- 
tra, use has been made of the so-called deep-water approximation to the linear wave dis- 
persion relationship. In general terms, the assumption is made that the water through 
which the waves are propagating is sufficiently deep to make the wave phase speed a func- 
tion of wave number alone, rather than a function of both wave number and water depth. 

Wave spectra to be used for refraction model verification, however, must necessarily be 
observed in nearshore regions in which the water depth is relatively shallow. For such 
applications, use of deep-water relationships for transforming apparent spectra observed 
from an aircraft platform will result in distorted estimates of the true nearshore spectra. 

Transformation equations valid for all depth zones are developed in the present 
paper under the assumption of a unidirectional wave field. The expressions are developed 
in terms of the linear wave dispersion relationship and the wave group speed. An iterative 
solution to the transformation equations is outlined for a specified value of average water 
depth in the area in which the apparent spectrum is observed. Sample cases are presented 
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to illustrate the effect of varying the assumed water depth in the transformation of refer ~ 
ence theoretical apparent spectra. 


SYMBOLS 


true wave phase speed, meters/second 


apparent wave phase speed observed from aircraft platform, meters/second 


local water depth, meters 


average water depth in area in which apparent spectrum is observed, meters 
incremental energy density within a narrow band of the spectrum, meters^ 
acceleration of gravity, meters/second^ 


wave number, meter“l 

L 


k* nondimensional wave number, kd 

L wavelength, meters 

S(cu) true spectral density, or wave energy density per unit true radian frequency, 

meters^ 
radian/second 

^a(^a) apparent spectral density, or wave energy density per unit apparent radian 

frequency, — meters 

radian/second 

Sa,r(^a>W) reference apparent spectral density for specified windspeed W, defined 

by equation (19), meters^ 

radian/second 

Spm(^;W) Pierson-Moskowitz (ref. 8) spectral density for specified windspeed W, 

defined by equation (18), ^^ters j — ^ 

radian/second 

V aircraft platform speed, meters/second 
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W windspeed, meters/second 

^ positive root of quadratic equation (eq. (A7)) defining deep-water approximate 

frequency transformation, radians/second 

I root of transcendental frequency transformation equation (eq. (6)), meters"^ 

o) true radian wave frequency, radians/second 

w* nondimensional true wave frequency, oj 

ct)a apparent radian wave frequency, radians/second 

TRANSFORMATION OF APPARENT SPECTRUM 

Developed in this section are equations which relate the apparent wave spectrum 
observed from an aircraft sensor platform to the spectrum which would result from 
observation of the same wave field from a stationary platform such as a wave staff. It is 
assumed that a digitized record of sea-surface elevation measurements made by the air- 
borne sensor is treated by standard time -series analysis methods to produce a spectrum 
whose ordinate values, when multiplied by (g/2)(Water density), yield an estimate of the 
distribution of wave energy density over a range of positive apparent frequencies, or fre- 
quencies observed from the moving aircraft platform. Such a treatment of the data record 
inherently assumes that the length of the record is sufficiently short and, thus, the varia- 
tion in water depth along the aircraft flight path is sufficiently small, to ensure homoge- 
neity of the data sample. For simplicity it is also assumed that the observed wave field 
is unidirectional (i.e,, all waves are traveling in the same direction) and that the aircraft 
is flying in that same direction at a constant speed which is greater than that of the fastest 
wave observed. 



Apparent Frequency Transformation 

Under the assumption of a unidirectional wave field being observed from an aircraft 
flying in the direction of wave propagation, the apparent phase speed of a given wave, or 
the phase speed relative to the airborne sensor, is 

Ca = c - V (1) 


4 



where Ca is the apparent phase speed, c is the true phase speed, or phase speed seen 
from a stationary platform, and V is the aircraft speed. The phase speed can be defined 
in terms of radian wave frequency and wave number by the equation 



where cu is the radian wave frequency and k is the wave number, equal to 27 t/L. 

Since wavelength (and, thus, wave number) is the same in both the stationary and the mov- 
ing coordinate systems (ref. 4, p. 302), equation (2) relates both true phase speed to true 
frequency and also apparent phase speed to apparent frequency. Equation (1) can be re- 
written then in terms of radian frequencies as 

== cu - kV (3) 

where u> is the true wave frequency and coa is the apparent wave frequency. Since the 
aircraft speed is assumed to be greater than the phase speed of any given wave being 
observed, the apparent frequency (see eq. (1)) is always negative. In order to relate the 
transformation to positive apparent frequencies which result from time -series analysis 
of the actual data record from the airborne sensor, it is convenient to redefine the fre- 
quency transformation as the negative of equation (3), that is, 

coa = kV - CO (4) 

Frequency and wave number are also uniquely related, according to linear wave 
theory, by the dispersion relationship 

co2 = gk tanh kd (5) 

where g is the acceleration of gravity and d is the local depth of the water through 
which the given wave is propagating. By assuming that the local depth in equation (5) can 
be replaced by an average depth over the observation area and by combining equations (4) 
and (5), the apparent frequency can be defined in terms of wave number as 

kV - (gk tanh kd)^/^ - coa = 0 (6) 

where d is the average water depth in the area in which the apparent spectrum is 
observed. Equations (5) and (6) in combination, then, relate a given apparent frequency 
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to its corresponding true frequency (or frequency observed from a stationary platform) 
for specified values of V and d. 

Spectral Density Transformation 

The apparent energy density spectrum derived from time -series analysis of the 
airborne measurements shows the distribution of energy density (energy per unit ocean 
surface area) with respect to apparent wave frequency. The incremental energy density 
AE within a narrow frequency band of the apparent spectrxun can be written as 

AE = Sa(coa|i) (waji^j - c^a|i) (7) 

where a>a|j and are the limits of the discrete apparent frequency band and 

Sa(<^a|^) is the apparent spectral density in this band. The corresponding band in the 
true spectrum is representative of the same physical waves, viewed from a stationary 
platform. Thus an equivalent amount of energy density must reside in the corresponding 
band in the true spectrum, that is, 

AE = S(co|i) - wli) (8) 

where w|i and are the transformed values of o)a|j and u>a|i+i> and S(a>|j) is 

the true spectral density in the transformed band. Then, setting equation (8) equal to 
equation (7) results in 

S(‘^li) (w|i+i - w|i) = Sa(c^a|i) (^a|i+i - ^a|i) (9) 

In the limit as the width of the frequency band approaches zero, equation (9) becomes 
S(o)) dw = Sa(wa) dcoa 


or 


S{oj) = Sa(wa)-^ 
dw 


The derivative dwa/dw can be determined directly from equation (4) as 


da>a _ y ^ _ 1 = V 
dw dw do)/dk 


(10) 


( 11 ) 


6 



Equation (10) can then be rewritten as 


S(w) = Sa(wa) 


' V 
da>/dk 



( 12 ) 


In physical terms, the derivative dw/dk is known as the wave group speed (ref. 7, p. 187), 
or the speed at which energy is transmitted by a wave train. Thus the energy density 
transformation is, by equation (12), a function of the ratio of aircraft speed to the speed 
at which energy is transmitted by the wave train being observed. 

A general expression for dcu/dk can be found by differentiating the dispersion 
relationship, equation (5), as 

■4^ = -^(tanh kd + kd sech^ kd) (13) 

dk 


so that the spectral density transformation can be written in final form as 

S{oo) = Sa(o^a) I _ A (14) 

\tanh kd + kd sech^ kd / 

where co, 0 ;^, and k are functionally related by equations (5) and (6) and d, again, is 
the average depth in the area in which the apparent spectrum is observed. 

Solution to Transformation Equations 

It was found under the assumptions stated in the previous sections that the transfor- 
mation of an apparent energy density spectrum, observed by an aircraft sensor, into the 
true energy density spectrum which would result from observation of the same wave field 
from a stationary platform can be accomplished by use of the following three equations: 

kV - (gk tanh kd)^' ^ = 0 (6) 


S(o;) = Sa(a)a) 


2coV/g 


— — 2 ~ 
tanh kd + kd sech kd 


- 


(14) 


= gk tanh kd 


(15) 
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Because of their transcendental nature, a closed-form solution to the transformation equa- 
tions cannot be found. The studies presented in references 4 to 6 circumvented this prob- 
lem by using what is known as the deep-water approximation. The assumption is made 
that the water depth is sufficiently large that for all wave numbers, tanh kd - 1 and 
sech kd 0, With this deep-water assumption the three transformation equations can be 
reduced to closed-form expressions which are given in the appendix. For application to 
apparent spectra observed in waters of relatively shallow depth, however, use of these 
deep-water approximations can result in substantial error. The magnitude of this error, 
along with the error which could result from certain other approximate solutions to the 
transformation equations, is discussed in the appendix. 

It is advisable, then, to seek a general solution to the transformation equations valid 
for all depth zones, A general solution can be found by first considering equation (6), the 
apparent frequency transformation in terms of wave number. For specified values of V, 
and d, standard iteration techniques can be used to find a root k = ^ of equation (6) 
which is unique under the assumption that the aircraft speed is greater than that of the 
fastest wave observed. The true frequency corresponding to the specified apparent fre- 
quency can then be found by using equation (15) as 


= (gl tanh 


and the true spectral density can be found by using equation (14) as 


s[co(d] = Sa(wa) 


2 co(g) V/g 

tanh + |d sech^ |d 


(16) 


(17) 


RESULTS AND DISCUSSION 


Sample calculations have been made which illustrate the effect on the transformation 
of apparent spectra of varying the average depth in the area in which the apparent spectra 
were observed. The theoretical fully developed wave spectrum form proposed by Pierson 
and Moskowitz (ref. 8) was used to form a reference apparent spectrum as follows; 

(1) The unidirectional wave spectrum as a function of true frequency for a given 
windspeed was specified as the Pierson-Moskowitz spectrum, 


Spm(^)W) = 8.1 X 10 ^ 


2 - exp 


-0.74 



4, 


(18) 
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where W is the speed of the wind generating the fully developed wave field. 

(2) A closed-form expression for the reference apparent spectrum for windspeed 
W was derived by applying the deep-water approximate transformation equations given 
in the appendix in inverse fashion to equation (18). For a specified aircraft speed V, the 
reference apparent spectrum for windspeed W was then derived as 


Sa,r(^a;W) - 


SpM(<^;W) 

2coV/g - 1 


(19) 


where 


CO 



CO 


(20) 


as shown in the appendix. 

The reference apparent spectrum at windspeed W, given by equations (19) and (20), was 
then retransformed by using the general solution of the transformation equations (eqs. (6), 
(14), and (15)) for a specified value of average water depth in the observation area. 

Reference apparent spectra for an aircraft speed V = 75 m/sec and various wind- 
speeds (10, 15, and 20 m/sec) were transformed into true spectra for assumed average 
water depths of 25 m, 50 m, 100 m, and deep water (d = «^). Comparison of the trans- 
formed spectra for W = 10 m/sec, given in figure 1, indicates that the spectra for water 
depths of 50 m, 100 m, and deep water are nearly identical. In the 25 -m transformed 
spectrum, there is a smaller value of the maximum spectral density accompanied by a 
slight broadening of the prominent portion of the spectrum. 

True spectra obtained by transforming the reference apparent spectrum for a wind- 
speed W = 15 m/sec are compared in figure 2. The transformed spectra for deep water 
and d = 100 m are again nearly identical. The 50-m transformed spectrum shows a 
smaller value for the maximum spectral density and a slight broadening in the prominent 
portion of the spectrum. The 25-m transformed spectrum shows an even smaller value 
for the maximum spectral density, a more pronounced broadening in the prominent portion 
of the spectrum, and a downward shift in the frequency at which the maximum spectral 
density occurs. 

True spectra obtained by transforming the reference apparent spectrum for a wind- 
speed W = 20 m/sec are compared in figure 3. Of particular interest in figure 3 is the 
shift to a significantly lower frequency range of the prominent portion of the 25-m trans- 
formed spectrum compared with the location of the prominent portions of the other spec- 
tra. In addition, the maximum spectral density is greater in the 25-m spectrum than in 



the 50-m spectrum, a result contrary to the trend in figures 1 and 2. The disparity among 
the transformed spectra shown in figure 3 indicates that when the generating windspeed is 
high and, by inference, much of the energy in the wave field is expected, a priori, to reside 
at relatively low true frequencies, the use of a proper depth value in transforming an 
apparent spectrum is vitally important. 

It should be reemphasized that, for a particular windspeed, the area under a trans- 
formed spectrvim, or total energy density, is invariant with assumed water depth (as seen 
from eqs. (7) to (9)), Thus, while preserving the total energy density of the wave field, 
use of an incorrect depth value in the transformation of the apparent spectrum of interest 
can lead to quite distorted estimates of the true spectral distribution of energy density. 

CONC LU DING RE MARKS 

The present paper has considered the transformation of a unidirectional apparent 
ocean wave spectrum observed from an aircraft sensor platform into the corresponding 
true spectrum, or the spectrum which would result from observation of the same wave 
field from a stationary platform. General transformation equations were developed in 
terms of the linear wave dispersion relationship and the wave group speed. An iterative 
solution to the transformation equations was outlined for specified values of aircraft speed 
and average water depth in the area in which the apparent spectrum is observed. 

Results obtained by using the iterative solution to transform sample theoretical 
apparent spectra indicate that significantly different transformed spectra can arise from 
a single apparent spectrum transformed with different assumed values of average water 
depth. For experimental cases in which much of the energy in a wave field can be 
expected, a priori, to reside at relatively low frequencies, use of an incorrect value for 
water depth in the transformation of the apparent spectra of interest can lead to quite dis- 
torted estimates of the true spectral distribution of wave energy. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
May 4, 1976 


10 



APPENDIX 


APPROXIMATE SOLUTIONS TO TRANSFORMATION EQUATIONS 

Approximate solutions to the equations defining the transformation of a unidirec- 
tional apparent wave spectrum into its corresponding true spectrum are examined in this 
appendix. A discussion of the errors associated with these approximate solutions serves 
to emphasize the need for an iterative solution to the exact transformation equations. The 
exact equations valid for any depth zone are stated as 

kV - w - o)a = 0 (Al) 

- *) 

where co and k are related by the linear wave dispersion relationship 

= gk tanh kd (A3) 

and do)/dk, which in physical terms is the wave group speed, is found by differentiating 
equation (A3) as 


= #-(tanh kd + kd sech^ kd) (A4) 

dk 2w 

By assuming that the water in the region in which the apparent spectrum is observed 
is sufficiently deep, so that for all wave numbers, tanh kd ~ 1 and sech kd ~ 0, equa- 
tions (A3) and (A4) can be approximated by 

« gk (A5) 


and 


dk 2o) 


(A6) 


Substitution of k as defined by equation (A5) into equation (Al) results in a closed-form, 
quadratic equation for co in terms of Wa which can be called the deep-water approxi- 
mate frequency transformation: 
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APPENDIX 


- w - wa = 0 (A7) 

g 

The positive root w = ? of equation (A7) can be substituted for w into equation (A6), 
and insertion of that expression for dw/dk in equation (A2) results in the deep-water 
approximate spectral density transformation: 

S(?) - Sa(o.a) - ij (A8) 

Use of these approximate relationships in the transformation of apparent spectra 
observed in regions of relatively shallow water can, however, result in substantial error. 
This error can best be examined by considering nondime nsional versions of the disper- 
sion relationship and the wave group speed. By defining co* = d/g and k* = kd, the 
dispersion relationship and group speed given in exact form by equations (A3) and (A4), 
respectively, can be written in nondimensional form as 

(o)*)2 = k* tanh k* (A9) 


and 


du>* _ tanh k* + k* sech^ k* _ tanh k* + k* sech^ k* 

2(k* tanh k*)^/2 


(AlO) 


Similarly, the deep-water approximate versions of the dispersion relationship and group 
speed, given by equations (A5) and (A6), respectively, can be written in nondimensional 
form as 


(Ct)*)2 s: k* 


(All) 


and 


do)* ^ 1 ^ 1 

dk* 2o)* 2(k*)^/^ 


(A12) 


Plots of the exact and deep-water approximate versions of the nondimensional dispersion 
relationship (eqs. (A9) and (All)) as a function of k* are compared in figure 4, and sim- 
ilar plots of the nondimensional group speed (exact version given by eq. (AlO), and deep- 
water approximate version given by eq. (A12)) are compared in figure 5. It can be seen 
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in these figures that for small values of k* (less than about 2,0), there is considerable 
difference between the exact and approximate versions of both the dispersion relationship 
and the group speed. Thus, the use of the deep-water approximation in the transformation 
of an apparent -frequency — spectral -density pair when k* = kd is small will result in 
both a misplacement of the transformed pair with respect to true frequency and a distor- 
tion with respect to true spectral density. 

In consideration of the disagreement at small values of k* between the exact dis- 
persion relationship and its deep-water approximation, the temptation exists to proceed 
to a piecewise -continuous approximation in which distinct functions are defined to approx- 
imate the exact dispersion relationship in adjacent ranges of k*. Candidate functions for 
such an approximation can be deduced by inspection of the exact nondimensional dispersion 
relationship shown in figure 4. For small values of k* (k* < 0.5), w* is seen to vary 
in a nearly linear fashion with k*. This follows from equation (A9), since tanh (k*) = k* 
for small values of k* (shallow water), so that ~ (k*)^. For k* > 2.0, the deep- 

water approximation (w*)2 » k* can obviously serve as the approximating function. For 
values of k* between 0.5 and 2.0, a transition function of the form, for example, 

~ aQ + aj(k* - 0.5) 

where aQ and a^ are constants, could be chosen so as to intersect the so-called 
shallow -water approximation at k* = 0.5 and the deep-water approximation at k* = 2.0, 
Shown in figure 6 is this prototype piecewise -continuous approximation, which, stated in 
equation form, is 


|^(k*)^ for k* < 0.5 

(a)*)2 = <Jo.25 + 1.167(k* - 0.5) for 0.5 g k* s 2.0 (A13) 

k* for k* > 2.0 

The maximum error in frequency which would be incurred by using this approximation is 
less than 10 percent. However, there would be a more severe problem encountered with 
the use of the corresponding approximation to the group speed dco*/dk* which, from 
equation (A13), is 


dco* 

dk* 


0.5833 


[0.25 + 1.167(k* - 0.5)]^/^ 


2(k*) 


1/2 


for k* < 0.5 
for 0.5 g k* g 2.0 

for k* > 2.0 


(A14) 
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APPENDIX 


A plot of equation (A 14) is compared with a plot of the exact nondimens ional group speed 
(given by eq. (AlO)) in figure 7, There are substantial differences between the two curves 
in certain ranges of k*, in addition to discontinuities in the approximate group speed at 
the points of intersection of the individual functions in the piecewise -continuous approxi- 
mation to the dispersion relationship. Sharp changes would thus be induced in the trans- 
formed spectra through use of such piecewise -continuous approximations to obtain a 
closed-form solution to the apparent -spectrum transformation equations. 
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Figure 2.- Transformed spectra for windspeed of 15 m/sec and various average water 

depths. Aircraft speed V = 75 m/sec. 
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Figure 3.- Transformed spectra for windspeed of 20 m/sec and various average water 

depths. Aircraft speed V = 75 m/sec. 




.2 


L.l 



J I I 1 ! 1 I 

.5 1.0 1.5 2.0 2.5 3.0 3.f 

k* 


Figure 5.- Exact nondimensional wave group speed du)*/dk* and 
deep-water approximation. 







Figure 7.- Exact nondimensional wave group speed and expression derived from 
piecewise -continuous approximation to dispersion relationship. 
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